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Abstract We investigate strong-coupling corrections to the specific heat Cy in the normal 
state of an ultracold Fermi gas in the BCS-BEC crossover region. A recent experiment on 
a ®Li unitary Fermi gas [M. J. H. Ku, et. al. Science 335, 563 (2012)] shows that Cy is 
remarkably amplified near the superfluid phase transition temperature T^, being similar to 
the well-known d-structure observed in liquid ^He. Including pairing fluctuations within 
the framework of the strong-coupling theory developed by Nozieres and Schmitt-Rink, we 
show that strong pairing fluctuations are sufficient to explain the anomalous behavior of Cy 
observed in a ®Li unitary Fermi gas near T^. We also show that there is no contribution from 
stable preformed Cooper pairs to Cy at the unitarity. This indicates that the origin of the 
observed anomaly is fundamentally different from the case of liquid ^He, where stable ^He 
Bose atoms induce the d-structure in Cy near the superfluid instability. Instead, the origin 
is the suppression of the entropy S, near Tc, due to the increase of metastable preformed 
Cooper pairs. Our results indicate that the specific heat is a useful quantity to study the 
effects of pairing fluctuations on the thermodynamic properties of an ultracold Fermi gas in 
the BCS-BEC crossover region. 

PACS numbers: 03.75.Hh, 05.30.Fk, 67.85.Lm. 


1 Introduction 

Ultracold Fermi gases have recently gathered much theoretical and experimental attention as 
highly controllable many-body system s A tunable pairing interaction associated 

with a Feshbach resonanccii^ has realized the so-called BCS (Bardeen-Cooper-Schrieffer)- 
BEC (Bose Einstein condensation) crossover phenomenon, where we can study a Fermi 
superfluid and a Bose superfluid in a unified manner. In the BCS-BEC crossover region, 
the system is dominated by strong pairing fluctuations, so that the formation of preformed 
Cooper pairs, as well as their effects on physical properties, have extensively been dis- 
cussed>2ii2iii. Recently various thermodynamic quantities have become accessible in the 
field of cold Fermi gas physic s 12,13,14,15 ^ allowing us to directly compare theoretical stud- 


1 iDepartment of Science and Technology, Keio University, Japan 







2 


P. van Wyk* et al. 


ies of the strong-coupling properties of an ultracold Fermi with experimental 

results. 

Since the formation of preformed Cooper pairs due to a strong pairing interaction is 
expected to affect the entropy S of the system, the specific heat (which is directly related 
to the change in entropy) is a useful thermodynamic quantity to study this key many-body 
phenomenon in the BCS-BEC crossover region. Indeed, a recent experiment on a ®Li unitary 
Fermi gas— observed that the specific heat at constant volume Cy deviates from the linear 
temperature dependence known in a normal Fermi liquid—. In the normal state near Tc, 
the observed Cy exhibits an anomalous enhancement with decreasing temperature, which 
resembles the so-called d-structure observed In liquid ^HenSS. Although this resemblance 
infers that some form of Bose excitations are responsible for the anomalous behavior of Cy 
in a ®Li unitary Fermi gas near T^, it is unclear to what extent a unitary Fermi gas with 
strong pairing fluctuations can be regarded as an interacting Bose system, like liquid ^He. 

In this paper, we investigate the specific heat at constant volume Cy, and effects of 
pairing fluctuations in the BCS-BEC crossover regime of an ultracold Fermi gas, above the 
superfluid phase transition temperature Tc- Including pairing fluctuations within the frame¬ 
work developed by Nozieres and Schmitt-Rink (NSR)^, we show that strong-coupling cor¬ 
rections to Cy, near Tc, in the unitarity limit are sufficient to explain the anomalous en¬ 
hancement of Cy that has recently been observed in a ®Li Fermi gas. We also clarify that 
this phenomenon originates from pairing fluctuations leading to the formation of metastable 
preformed pairs, which is quite different from the case of liquid ^He, where stable ^He Bose 
atoms induce the d-structure in the temperature dependence of Cy near the superfluid insta¬ 
bility—. Throughout this paper we take h = k^ = \, and the system volume V is taken to be 
unity, for simplicity. 


2 Formulation 


We consider a two-component Fermi gas, described by the BCS Hamiltonian, 

H ='^^pcl,crCp,a--U ^ ‘^p+,/2,t‘^-p+g/2,i‘^-P'+?/2d‘^P'+¥/2,T> 

p<^ p,p',q 


( 1 ) 


where c^.o- is the annihilation operator of a Fermi atom with pseudospin cr =1, J,, describing 
two atomic hyperfine states, = Ep-I{2m )—is the kinetic energy, measured from 
the Fermi chemical potential p, where m is an atomic mass. The pairing interaction -U (< 0) 
is assumed to be tunable by adjusting the threshold energy of a Feshbach resonance^. As 
usual, we measure the interaction strength in terms of the inverse j-wave scattering length 
a^*, which is related to -U through 


AnUs —U 

m 1 _ {/ y 

i ^ djp ISp 


(2) 


We include fluctuations in the Cooper channel within the framework of the strong¬ 
coupling theory developed by Nozieres and Schmitt-Rink (NSR)^’. In the Green’s function 
formalism, the NSR theory is described by the single-particle thermal Green’s function. 


Gp{iw„) = + G“(iw„)Z-p(iw„)G°(/w„), 


( 3 ) 
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Fig. 1 Self-energy Zg{iojn) in the NSR Green’s function in Eq. s. The solid line and dashed line are the 
non-interacting Green’s function G®, and the pairing interaction —U, respectively. 


where ai,, is the fermion Matsubara frequency, and = iai„ —^p is the Green’s func¬ 

tion for a free Fermi particle. The NSR self-energy Sp{i<jj„), which describes strong-coupling 
corrections to single-particle excitations, is diagrammatically drawn as Fig.[T] which gives, 

Ep(icj„) = r ^ rg(iv„)Gp_^(iLiJn - iv„). (4) 

qJvn 


Here, is the boson Matsubara frequency. The NSR particle-particle scattering matrix, 

1 




Anas 


riqiiVn) 2 e„ 


describes fluctuations in the Cooper channel, where 

1 -/(W2)-/(^P-.7/2) 

~^p+ql2 ~^p-q/2 


nq{lV„) = - > -^- 


( 5 ) 


( 6 ) 


is the lowest order pair-correlation function (where f{x) is the Fermi distribution function). 
We briefly note that, although the pair correlation function Ilqiiv,,) has an ultraviolet diver¬ 
gence, this singularity is renormalized into the j-wave scattering length in Eq. l|2]i. 

The specific heat at constant volume Cy is calculated from the formula. 



V,N 


Here, the NSR internal energy E is given by 


( 7 ) 


E = 2tJ^ 

pMi 


^p "f ^22piiojfi) 


G p(iCi>n ) 


— T Eq(iv„) 


T -^nqiiVn) +H^nq{iVn) 
di 0)1 


( 8 ) 


where we have retained only terms to 0(Zp{ia>„)), so as to be consistent with the NSR 
Green’s function in Eq. ([3ll. In evaluating l|7]i and ([8]l, we need the value of the chemical 
potential which is obtained from the equation for the number N of Fermi atoms, 

N = 2r y Gp(m„). 

p,iLO„ 


( 9 ) 
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Fig. 2 (Color online) (a) Compaiison of our result with the recent experiment on a ®Li unitary Fermi gasJ^. 
Cy shows the specific in a free Fermi gas. (b) Calculated specific heat Cv in the BCS-BEC crossover regime 
of an ulteacold Fermi gas above 


Equation (O is also used in self-consistently determining the superfluid phase transition 
temperature Tc from the Tc -equation— , 



which is obtained from the Thouless criterion—, which states that the superfluid instability 
occurs when the particle-particle scattering matrix /’g(iV„) has a pole at ^ = 0 and v„ = 0. 

In our numerical calculations, we tune the interaction strength by adjusting the value of 
the dimensionless parameter . We briefly note that, in the unitarity limit ((A:Faj)“* = 

0), the interaction effect formally disappears in Eqs. © and ( fTOll . in the sense that the factor 
mHAnas) in these equations vanishes. 


3 Specific heat in the BCS-BEC crossover region above Jc 

Figure |2ja) shows the specific heat Cy in a unitary Fermi gas above Tq. Starting from the 
high temperature region, Cy decreases from the classical result Cy = {3I2)N, with decreas¬ 
ing temperature. Since this high temperature behavior is well described by the case of a 
free Fermi gas (Cy), this temperature dependence may be simply attributed to the Fermi 
statistical effect. While the free Fermi gas result Cy monotonically decreases with decreas¬ 
ing temperature, the specific heat Cy at the unitarity is enhanced in the low temperature 
region, Tc <T <2Tc. As seen in Fig.j^a), this enhancement semi-quantitatively agrees with 
the recent experiment on a ®Fi Fermi gasi^iS. Since the physical properties of a unitary 
Fermi gas are known to be dominated by strong pairing fluctuations near 8^9^10,16^17,21^24 ^ 
this anomaly in Cy is also considered as a strong-coupling phenomenon associated with 
strong fluctuations in the Cooper channel. Indeed, as shown in Fig. |2jb), this enhancement 
is remarkable in the unitary regime. 
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At a glance, the anomalous enhancement of the specific heat Cv seen in Fig.|2a) resem¬ 
bles the case of liquid ^Hci^S, where stable atoms constituting this Bose liquid are respon¬ 
sible for the famous /i-structure in the temperature dependence of the specific heat near the 
superfluid instability. Since pairing fluctuations in a unitary Fermi gas are accompanied by 
preformed Cooper pairs, it is an interesting problem whether or not the observed enhance¬ 
ment of Cv can be understood by regarding a unitary Fermi gas as a “Bose gas” of stable 
preformed Cooper pairs. 

To examine this, it is convenient to divide the number equation l|9]l into the sum of the 
free fermion part, 

Af = 27’ E G°(im„) = 2j]/(^p), (11) 

p,iCL)n p 

and the fluctuation contribution coming from the second term in Eq. j[3]l, 

Afl = -7’ E rqiiVn)-^nq{iVn) = C dcUHe(m)pB(w). (12) 


Here, nB(m) is the Bose distribution function, and 


Pb(oj) = y — Im 
¥ 


Ca(!V„ ^ <jj + i6)—n„(iv„ -> m-fid) 
op 


(13) 


can be interpreted as the “molecular density of states” (where 6 is an infinitesimally small 
positive number). Equation lll2b implies that Afl is related to the bosonic character of the 
system. Evaluating the contribution from real poles (Qq’s) of the analytic-continued particle- 
particle scattering matrix /’^(;V„ —> ai + id) to Eq. (112b (= Nst), one ha s^'i^^ 


a^st= y 

Qq 


g^^nqiQq) 

■^IJqiQq) 


nsi^q). 


(14) 


where the summation is taken over real poles (Qq’s). Apart from the prefactor of the Bose 
distribution function in Eq. ( 114b . Ast is found to be directly related to the number of stable 
molecular bosons having the dispersion Qq. In the extreme BEC limit (where the system is 
well described by a gas of tightly bound molecules), the prefactor is reduced to two, and 
the molecular dispersion becomes Qq = jlM (where M = 2m), as expected. To conclude, 
the fluctuation contribution Afl = iVsT + Nsc can be decomposed into twice the number of 
stable preformed Cooper pair molecules Asx and the so-called scattering part Asc ^*'^^ ; the 
latter of which physically describes the number of metastable preformed pairs with a finite 
lifetime. 

As shown in Fig. [3] while the fluctuation contribution Afl at is dominated by stable 
molecules (Asx) in the strong-coupling BEC regime, they are absent in the unitarity limit^. 
This indicates that, although the enhanced specific heat Cv(T ~ Tc) shown in Fig. |2la) is 
similar to the d-structure in the case of liquid ^He, their origins are different from each 
other. While stable '^He atoms are responsible for this anomaly in a liquid ^He, the sup¬ 
pression of the entropy due to the formation of metastable preformed pairs, with decreasing 
temperature, leads to the amplification of the specific heat near Tc in a unitary Fermi gas, 
through the thermodynamic relation. 


Cv = T 




V,N 


( 15 ) 
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Fig. 3 (Color online) Calculated number of stable preformed pairs Nst(T = Tc), as well as and contribution 
of scattering states N^ciT = Tc). The fluctuation contribution A^fl is given by the sum of the two as Wfl = 
Nsc + iVsT—■ 


We see in Fig.[2]that the number Ngc of metastable Cooper pairs at first increases, as 
one approaches the unitarity regime from the weak-coupling BCS side. Since the increase 
of this fluctuation contribution suppresses the entropy S{T ~ Tc), the specific heat Cy is en¬ 
hanced near Tc, as shown in Fig.|2b). Such a fluctuation effect on the entropy S is, however, 
suppressed in the strong-coupling BEC regime, because metastable preformed Cooper pairs 
are replaced by stable tightly bound molecules with a finite binding energy— Ebind = llma^. 
(See Fig. [3]) This is the reason for the non-monotonic behaviour of Cv{T = Tc) with respect 
to the interaction strength seen in Fig.|2b). 

In the strong-coupling BEC regime, the system is well described by a gas of A^b = N/2 
tightly bound molecules with the molecular mass M = 2m. As a result, Cy increases with 
decreasing the temperature near Tc when 2"0.8, as in the case of an ideal Bose gas 

(although we do not explicitly show this in this paper). In the BEC limit, Cv{T = Tc) just 
equals the specific heat =; 1.93Ab of an ideal molecular Bose gas at the BEC phase 
transition temperaturei^^. 


4 Summary 

To summarise, motivated by the prospective usefulness of the specific heat Cy as a probe to 
study effects of strong pairing fluctuations in an ultracold Fermi gas, we have discussed 
strong coupling corrections to this thermodynamic quantity in the BCS-BEC crossover 
region. Within the framework of the strong-coupling theory developed by Nozieres and 
Schmitt-Rink, we showed that Cy is enhanced in the unitary regime near Tc. We also showed 
that this anomalous enhancement is sufficient to explain the recent experiment on a ®Li uni¬ 
tary Fermi gas, where the observed specific heat is amplified with decreasing the temperature 
near Tc- 

Although this behaviour looks similar to the T-structure of the temperature dependence 
of the specific heat in a liquid ^He near the superfluid phase transition temperature, we 
pointed out that the origin of the former is fundamentally different from the latter. Indeed, we 
found that there are no stable preformed Cooper pairs at the unitarity, which is quite different 
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from the latter case, where “^He atoms are always stable. Instead, the bosonic character of 
a unitary Fermi gas near Tc is dominated by the increase of metastable preformed pairs as 
the temperature is lowered. This naturally suppresses the entropy near Tc, leading to the 
amplification of Cy. Our result indicates that the anomalous amplification of the specific 
heat observed in a ®Li unitary Fermi gas can be fully explained by fluctuating Fermi pairs, 
without the inclusion of stable bosons. 

In this paper, we have only treated the normal state. Since the specific heat has also been 
observed in the superfluid phase, it is important to extend the present theory to the region 
below Tc. In addition, since a real Fermi gas is trapped in a harmonic potential, effects 
of spatial inhomogeneity on the specific heat is also a crucial future problem. Since the 
specific heat is a fundamental thermodynamic quantity, our results would contribute to the 
further understanding of thermodynamic properties of an ultracold Fermi gas in the BCS- 
BEC crossover region. 
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